Linear scaling between momentum- and spin scattering in graphene 
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Spin transport in graphene carries the potential of a long spin diffusion length at room tempera- 
ture. However, extrinsic relaxation processes limit the current experimental values to 1-2 ^m. We 
present Hanle spin precession measurements in gated lateral spin valve devices in the low to high 
(up to 10^^ cm~^) carrier density range of graphene. A linear scaling between the spin diffusion 
length and the diffusion coefficient is observed. We measure nearly identical spin- and charge dif- 
fusion coefficients indicating that electron-electron interactions are relatively weak and transport is 
limited by impurity potential scattering. When extrapolated to the maximum carrier mobilities of 
2x10^ cm^/Vs, our results predict that a considerable increase in the spin diffusion length should 
be possible. 

PACS numbers: 72.25.Hg, 73.63.-b 
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The high charge carrier mobilityJl|,i,[3 and spin dif- 
fusion length of micrometers [1, [^^measured at room 
temperature make graphene a possible candidate for 
future electronic and spintronic devices. This two- 
dimensional crystalline material has two electronic con- 
duction regimes, metallic where charge carriers are of 
one type (holes or electrons), and the region around the 
Dirac neutrality point where transport of electric current 
happens through small regions charged locally with holes 
or electrons ("electron- hole puddles"). The presence of 
such puddles yielding a finite local density |n| « 10^^ 
cm~^ was shown experimentally using scanning single 
electron transistor technique!^ and scanning electron 
spectroscopy 7, 8^ , with the intensity of the fluctuations 
being strongly enhanced by substrate impurities 

Experiments done so far on spin transport reveal room 
temperature spin relaxation times of the order of 100 
to 200 ps[|. This is well below the theoretically pre- 
dicted, intrinsic limit [lol| b ut might be explained if 
we consider extrinsic effects[ll|. Since hyperfine interac- 
tions at 300K are weak in graphitic systems, there are 
two possible mechanisms that can be held responsible 
for such a strong spin relaxationp^. scaling differently 
on the momentum relaxation. In case of the Elliot- Yafet 
mechanism (spin flip occurs with a finite probability at 
each momentum scattering center) the spin scattering 
time Ts is proportional to r^, the momentum scattering 
time, while the D'yakonov-Perel mechanism (spins pre- 
cess under the influence of local spin-orbit fields in be- 
tween scattering events) is characterized by cx t^^. To 
identify the scattering mechanism and find the ultimate 
limit on spin relaxation, one can thus investigate the link 
between spin transport and the electronic quality of the 
graphene, in particular the charge carrier mobility. Since 
the mobility is ill-defined at or in the vicinity of the Dirac 
neutrality point, we will use the diffusion coefficient de- 
fined as D = ^vpl and link it to the spin diffusion length 
Xs ~ \J Dts . Here vp is the Fermi velocity and I repre- 



sents the scattering mean free path. 

One way to study charge against spin diffusion would 
be comparing the results in a set of devices that dis- 
play significantly different carrier mobilities. However, 
it is experimentally challenging to fabricate consistently 
good ferromagnetic contacts to the graphene for such a 
set of samples. The option we choose here is to do the 
experiments on individual devices tuning the carrier den- 
sity from the metallic regime down to the lowest values 
and comparing the behavior of the spin transport to the 
changes in the charge diffusion coefficient. 

In this Rapid Communication we present a system- 
atic study of the spin transport and scattering at room 
temperature in single layer graphene samples on Si02 
substrate. The measurements are done at a wide range 
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1: (Color online.) a) SEM image of a spin valve de- 
vice. A 0.3 X 13fj,m strip with a cross shape in the middle 
(brown/light gray) was etched with oxygen plasma out of the 
original graphene flake (dark gray). The Co electrodes 1-10 of 
widths 90 to 800 nm and spacings 1 to 3.1 nm are also visible, 
b) Hanle precession measurements (dots) and flts (solid lines) 
at the Dirac neutrality point and in the metallic regime, with 
the injector/detector magnetization aligned parallel(P) and 
antiparallel( AP) . 
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of carrier densities with an accent on the Dirac neutral- 
ity point where the transport is difficult to model and 
Coulomb electron-electron interactions are expected to 
be the strongest 'l3'|. We compare this directly to the 
charge transport in the same samples to learn more about 
the diffusion phenomena and the interactions that lead 
to spin relaxation. 

The charge carrier transport in graphene in the metal- 
lic regime (at an energy E sufficiently far away from the 
Dirac neutrality point) can be described by a diffusion 
process characterized by the 2-dimensional charge diffu- 
sion coefficient D. The density of states (DOS) is given 



Vg(V) 

-30 -20 -10 10 20 30 40 50 



g.t,gs2TT\E\ 



(1) 



2) 

,-1 



with the twofold valley {gy = 2) and spin (gg 
degeneracies and the Fermi velocity vp ~ fO^ms 
By integration we can obtain the density n{Ep) = 
gvgsT^Ep/ {h'^Vp) with Ep the Fermi energy, and the Ein- 
stein relation a — e^vD allows for calculating the charge 
diffusion coefficient 
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Here e is the electron charge and Rg is the square resis- 
tance of the graphene layer, inverse of the conductivity 
(7. Finally, using the semiclassical Drude formula one 
can calculate the carrier mobility /i = {Rsne)~^ for the 
metallic regime. 

In order to determine the charge and spin diffusion 
coefficients experimentally we have fabricated field effect 
devices where the single layer graphene flake is contacted 
by several ferromagnetic cobalt electrodes. A scanning 
electron microscope (SEM) image of such a device is 
shown in Fig. [T^) . The graphene flakes are obtained by 
mechanical exfoliation from commercially available Kish 
graphite and deposited on a thermally oxidized, n+ + 
doped Si substrate (300 nm thick oxide layer). The Si 
substrate contacted by a bottom An electrode is used 
as electrostatic gate; applying a voltage Vg of typically 
tens of volts on it allows reaching carrier densities in 
the graphene \n\ < 10^^ cm~^ (electrons or holes) calcu- 
lated from the capacitance [iTj. A set of predefined Ti/Au 
markers help to accurately locate selected graphene flakes 
through optical- and atomic force microscope. The flakes 
are then etched with oxygen plasma into a cross shape, 
to allow for precise Hall type measurements using side 
contacts e.g. the ones labeled 5 and 6 in Fig. [IJi) as cur- 
rent injectors and contacts 4 and 7 as Hall voltage probes. 
The electrical contacts are patterned using electron beam 
lithography and evaporated thermally at a base pressure 
of ~ 10~^ mbar followed by standard lift-off technique. 
To achieve a high spin injection efflciencvp^ a 0.8 nm 
thick AI2O3 insulating layer was introduced between the 
graphene and the ferromagnet, resulting in contact resis- 
tances of the order 20-40 kil. 
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FIG. 2: (Color online.) Charge- and spin transport parame- 
ters plotted against charge carrier density and gate voltage, a) 
Square resistance; b) Non-local spin valve signal determined 
from spin valve- and from Hanle precession measurements; 
c) Charge- and spin diffusion coefficients; d) Spin relaxation 
time; e) Spin relaxation length. 



The graphene's square resistance i?s is determined 
from four-probe local measurements. Sending an elec- 
tric current from e.g. contact I to 4 and measuring the 
voltage drop between 2 and 3 we are sensitive only to 
the resistance of the graphene between contacts 2 and 3. 
Measuring the resistance against the applied gate volt- 
age (i.e. in function of carrier density) and normalizing 
it to the graphene length to width ratio yields the Rs 
curve plotted in Fig. [5] a). From such measurements we 
calculate the charge diffusion coefficient Z?c in the metal- 
lic regime using relation [2l see the solid line in Fig. [2j;). 
The decrease in carrier density comes with a decrease 
in the diffusion coefficient; the singularities in the calcu- 
lated Dc at the charge neutrality point will be discussed 
later. The asymmetry in the electron- versus hole diffu- 
sion at high densities visible in panel c) probably origi- 
nates from nonuniformities in the carrier density and can 
be traced back to the Rg measurement in panel a). Mea- 
surements of the Hall coefficient Rh (not shown) against 
the gate voltage using the cross contact geometry indi- 
cate the onset of the metallic regime at a carrier density 
n ~ ±0.5 X 10^^ cm~^ by displaying a clear f/n depen- 
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dence. The density value extracted from the Hall mea- 
surements in the metallic regime confirms the number 
calculated from the square resistance measurements and 
gate capacitance. 

The spin transport measurements are performed in 
the non-local geometry a spin-polarized current is in- 
jected e.g. through electrode 2 and extracted through 
electrode 1, while we measure the voltage between elec- 
trodes 3 and 4. There is no charge current flowing be- 
tween 3 and 4; the detected non-local signal Rni in an in- 
plane magnetic field is purely due to the effect of spins dif- 
fusing from the injector electrodes to the detectors. Sub- 
tracting Rni at parallel and antiparallel magnetic orien- 
tation of the injector/detector electrodes while scanning 
the gate voltage gives the spin valve signal Rsv that has 
a significant dependence on the charge carrier density as 
plotted in Fig.[2]b), solid line. The electrodes 1 and 4 are 
far enough not to contribute significantly to Rsvj there- 
fore we define L — 3.1/ini as the distance between the 
two inner electrodes. 

Applying a magnetic field Bz orthogonal to the sam- 
ple plane will result in Hanle spin precession. Measuring 
Rni while we sweep the magnetic field (i.e. we change the 
precession frequency) yields the curves in fig. [T]d) . Here 
two measurements are plotted for the metallic regime, 
Vg = +iOV (n ~ -2 X 10^2 cm-2), and two for the 
Dirac neutrality point, Vg = +9V^ with the central injec- 
tor/detector electrodes oriented parallel and antiparallel, 
respectively. The parallel-antiparallel signal difference at 
zero field is the same as the spin valve signal defined 
above and is plotted for different densities in Fig. Oa), 
dots. The advantage of a spin precession measurement 
is that it allows extracting the spin diffusion coefficient 
Ds and spin scattering time by fitting the measure- 
ments with the solutions to the I31och cquation[T2| for 
spin accumulation 
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where the first term on the left hand side describes the 
spin diffusion, the second term the spin relaxation and 
the last one the precession, with an effective Lande factor 
g = 2 and the Bohr magneton fiB- 

A set of precession measurements was done for differ- 
ent charge carrier densities; the resulting spin transport 
parameters Dg and are plotted in Fig.O:) and d). The 
spin diffusion length As = \/D^t^ is shown in panel e) . 

Examining Fig. 2, we see that D^, and all de- 
crease approximately by a factor of 2 when we approach 
the neutrality point. This results in a strong decrease in 
the detected spin valve signal as seen in Fig. [2)d) , consis- 
tent with the prediction of the formula for a four-terminal 
non-local spin injection geometry [l^: 
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where Wg = 300nm is the width of the graphene flake. 
The spin polarization of the injected current determined 
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FIG. 3: (Color online.) Spin- vs. cfiarge diffusion coefficient 
with tlie unbroadened DOS from Eq. [1] and the broadened 
version from Eq. [5] using a Gaussian broadening of FWHM~ 
176 meV. 



from this relation is P ~ 9%. Note that this value can be 
considered constant through the range of carrier densities 
we used, since the contact resistances span from 20 to 
40 kfl where impedance mismatch is suppressed [22|. 

Let us focus now on the diffusion of charge versus spin. 
As visible in Fig. [2t), in the high density case the values 
are practically identical for spin and charge. This is a 
striking observation, since the two physical entities, Dc 
and Ds, are determined from completely different types 
of experiments. 

However, for \n\ < 0.5 x lO^^cm^^ formula [T] yields un- 
physical values for the diffusion coefficient and results in 
a singularity at the Dirac neutrality point. This comes 
from the unrealistic assumption of vanishing carrier den- 
sity and DOS. To correct for it, one has to account for 
a broadened density of states v*{E) due to finite tem- 
perature, electron-hole puddles and possibly to the finite 
lifetime of electronic states. The simplest way to include 
all broadening effects in the DOS is to add a Gaussian 
broadening energy a in the form of 
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After integration we obtain 
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where erf is the Gaussian error function and the only 
undetermined parameter is the value of a. Replacing 
the DOS with the broadened version in formula [2] we 
plot the modified diffusion constant D* in function of 
the density together with the unmodified charge- and 
spin diffusion constants, see Fig. [3l We find good cor- 
respondence between D* and Ds both at low and high 
densities if and only if we choose an energy broadening 
of cr ~ 75meV [l^l, i.e. a Gaussian with full width at 
half maximum (FWHM) = 2^j2ln2a ~ 176 meV cor- 
responding to a density variation of An ~ ±0.7 x 10^^ 
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FIG. 4: (Color online.) Linear relationship between the spin 
relaxation length and the spin difltusion coefficient, extracted 
from Fig. 2 panels c) and e). 

cm^^. This is consistent with the literature values [7] 
attributed to electron-hole puddles in graphene on Si02, 
considering that our samples show a carrier mobility of 
only 3000 cm V Vs. [3 

The observation that the spin diffusion coefficient 
shows no considerable difference from the charge diffusion 
coefRcient indicates a minor role of Coulomb electron- 
electron interactions [lOjl. This is in agreement with the 
recent results of Li et al. j20l| where electron-electron inter- 
actions are detected in high carrier mobility, suspended 
graphene flakes only. The major mechanism for limiting 
the (spin) diffusion seems to be the impurity potential 



scattering. 

The most convincing argument comes, however, from 
the scaling between (spin- or charge) diffusion coeffi- 
cient and spin diffusion length. Plotting the values of 
against Dg as extracted from Fig. [2] shows a clear linear 
dependence for both the electron and the hole conduc- 
tion regime, see Fig. S) Since As = \JDts, the linear de- 
pendence means that the spin scattering time is directly 
proportional to the diffusion coefficient, i.e. to the mo- 
mentum scattering time. The experiments confirm thus 
an Elliot- Yafet type spin relaxation mechanism, in agree- 
ment with our earlier spin relaxation anisotropy studies 
presented in Ref. [2l|. 

In conclusion, we expect that improving the electronic 
characteristics of the graphene flake by e.g. removing 
the substrate (suspended graphene) , annealing with high 
electric currents and/or using selected starting material 
(higher purity graphite) shall both enhance the charge 
transport and prolong the spin scattering time. Assum- 
ing the Elliot- Yafet mechanism is still dominating at high 
carrier mobilities we can extrapolate the behavior shown 
in Fig. [4] to samples displaying charge carrier mobilities 
in the range of 2 x 10^ cm^/Vs reported recently, to pre- 
dict a possible room temperature spin diffusion length up 
to 100 micrometers [24]. 

This work is part of the research programme of the 
Foundation for Fundamental Research on Matter (FOM), 
which is financially supported by the Netherlands Or- 
ganisation for Scientific Research (NWO). We acknowl- 
edge the financial support of the Zernike Institute for 
Advanced Materials and NanoNed. 



[1] K. Bolotin, K. Sikes, Z. Jiang, M. Klima, G. Fudenberg, 

J. Hone, P. Kim and H. Stormer, Solid State Commun. 

146, 351 (2008) 
[2] S.V. Morozov, K.S. Novoselov, M.I. Katsnelson, F. 

Schedin, D.C. Elias, J. A. Jaszczak and A.K. Geim, Phys. 

Rev. Lett. 100, 016602 (2008) 
[3] X. Du, I. Skachko, A. Barker and E. Y. Andrei, Nat. 

Nanotechnol. 3, 491 (2008) 
[4] N. Tombros, C. Jozsa, M. Popinciuc, H.T. Jonkman and 

B.J. van Wees, Nature (London) 448, 571 (2007). 
[5] W. H. Wang, K. Pi, Y. Li, Y. F. Chiang, P. Wei, J. Shi, 

and R. K. Kawakami, Phys. Rev. B 77, 020402(R) (2008) 
[6] J. Martin, N. Akerman, G. Ulbricht, T. Lohmann, J. H. 

Smet, K. von Klitzing and A. Yacoby, Nat. Phys. 4, 144 

(2008) 

[7] A. Deshpande, W. Bao, F. Miao, C. N. Lau and B. J. 

LeRoy, Phys. Rev. B 79, 205411(2009) 
[8] Y. Zhang, V. W. Brar, C. Girit, A. Zettl and M. F. Crom- 

mie, Nat. Phys., to he published; arXiv:0902.4793 
[9] D. Huertas-Hernando, F. Guinea and A. Brataas, Phys. 

Rev. Lett. 103, 146801 (2009) 
[10] C.L. Kane and E.J. Mele, Phys. Rev. Lett. 95, 226801 

(2005) ; H. Min et al, Phys. Rev. B 74, 165310 (2006); 
D. Huertas-Hernando et al, Phys. Rev. B 74, 155426 

(2006) ; Y. Yao et al, Phys. Rev. B 75, 041401(R) (2007); 



M. Gmitra et aL, 'arXiv:0904.33J_5J (2009) 
[11] A. H. Castro Neto and F. Guinea, Phys. Rev. Lett. 103, 
026804 (2009) 

[12] J. Fabian, A. Matos-Abiague, C. Ertler, P. Stano and I. 

Zutic, Acta Phys. Slov. 57, 565 (2007); arXiv:0711.1461| 
[13] Y. Barlas, T. Pereg-Barnea, M. Polini, R. Asgari and A. 

H. MacDonald, Phys. Rev. Lett. 98, 236601 (2007); E. 

H. Hwang, Ben Yu-Kuang Hu and S. Das Sarma, Phys. 

Rev. Lett. 99, 226801 (2007) 
[14] M. Popinciuc, C. Jozsa, P.J. Zomer, N. Tombros, A. 

Veligura, H.T. Jonkman and B.J. van Wees, Phys. Rev. 

B, to be published 
[15] T. Ando, J. Phys. Soc. Jpn. 75, 074716 (2006) 
[16] F.J. Jedema et al, Nature (London) 416, 713 (2002) 
[17] A. K. Geim and K. S. Novoselov, Nature Mater. 6, 183 

(2007). 

[18] C. Jozsa, M. Popinciuc, N. Tombros, H.T. Jonkman and 
B.J. van Wees, Phys. Rev. Lett. 100, 236603 (2008) 

[19] CP Weber, N. Gedik, J.E. Moore, J. Orenstein, J. 
Stephens and D.D. Awschalom, Nature )London) 437, 
1330 (2005) 

[20] G. Li, A. Luican and E.Y. Andrei, Phys. Rev. Lett. 102, 
176804 (2009) 

[21] N. Tombros et al., Phys. Rev. Lett. 101, 046601 (2008) 
[22] We have considered the spin injection model for transpar- 



5 



ent contacts from [l^ that corrects formula |4] addressing 
the impedance mismatch. The influence of the correction 
on the injection efficiency is below 10%. 
[23] The value of a is chosen to cure the singularity at n « 
by fitting the broadened D* to Ds for densities |7i| < 
0.5 X 10^^ cm~^. At |n| > 3 X 10^^ cm"^ the broadened 
curve approaches the unbroadened one determined from 



the measurements. 
[24] In metals it is known that the ratio between spin- and mo- 
mentum relaxations is approximately the same for scat- 
tering induced by impurities or electron-phonon interac- 
tions. We assume that this is also the case in graphene. 



